Primordial large-scale electromagnetic fields from Gravitoelectromagnetic Inflation 
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We investigate the origin and evolution of primordial electric and magnetic fields in the early 
universe, when the expansion is governed by a cosmological constant Aq. Using the gravitoelec- 
tromagnetic inflationary formalism with Aq = 0, we obtain the power of spectrums for large-scale 
magnetic flelds and the inflaton field fluctuations during inflation. A very important fact is that our 
formalism is naturally non-conformally invariant. 
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^ ■ I. INTRODUCTION 

, ^ ' The origin of the primordial magnetic fields has been subject of a great amount of research[l|. The existence, 
strength and structure of these fields in the intergalactic plane, within the Local Superclusted, has been scrutinized 
recently 0- Many spiral galaxies are endowed with coherent magnetic fields of fiG (micro Gauss) strength [1, 0, H, 
S S [1] J having approximately the same energy density as the Cosmic Microwave Background Radiation (CMBR) . In 
particular, the field strength of our galaxy is -B ~ 3 x 10~^ G, similar to that detected in high redshift galaxies [9| and 
dam ped Lyman alpha clouds [l^l- Limits imposed by the high isotropy of CMB photons, obtained from the COBE 
^ I datadH restrict the present day strength of magnetic fields on cosmological scales to 10~^ G. It is very mysterious 
that magnetic fields in clusters of galaxies [i.e., on scales ~ Mpc], to be coherent • There are compelling indications 
of existence of large-scale microgauss magnetic fields in galaxy clusters. This would indicate that the entire universe 
is magnetized. There are two possible classes of mechanisms to produce cosmic fields depending on when they are 
generated: astrophysical mechanisms acting during large-scale structure formation, and mechanisms acting in the 
' primordial universe. The origin of these magnetic fields is not well understood yet. The seeds of these fields could 
^ ' be in the early inflationary expansion of the universe, when these fields were originated. The existence of primordial 
magnetic fields would affect both, the temperature and polarization anisotropics of the cosmic microwave background. 
It also provides a plausible explanation for the possible disparity between observations and theoretical fits to the CMB 
. power-spectrum. The ACBARjl^ and CBl[lJj experiments indicate continued power up to I ~ 4000, but WMAP 
• ' data predicts a rapidly declining power spectrum in the large multipole range [I^- This discrepancy is difficult to 
^ account from a returning of cosmological parameters. Among other possible explanations, an cosmological magnetic 
' field generated during inflation provides a plausible mechanism to produce excess power at high multipoles. Therefore, 
I the study of its origin and evolution in this epoch should be very important to make predictions in cosmology [l^. 
• During inflation the extension of the causally connected regions grows as the scale factor and hence faster than in 
the decelerated phase. This solves the horizon problem. Furthermore, during inflation the contribution of the spatial 
curvature becomes very small. The way inflation solves the curvature problem is by producing a very tiny spatial 
$H ' curvature at the onset of the radiation epoch taking place right after inflation. The spatial curvature can well grow 
during the decelerated phase of expansion but it will be always subleading provided inflation lasted for sufhciently long 
time. It is natural to look for the possibility of generating such a large-scale magnetic field during inflation. However, 
the FRW universe is conformal flat and the Maxwell theory is conformal invariant, so that magnetic field generated at 
inflation would come vanishingly small. Therefore, the conformal invariance must be broken to generate non-trivial 
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magnetic fields. Various conformal symmetry breaking mechanisms have been proposed so far 17|. Magnetogenesis 
has been studied also during the electroweak phase transition[18j. Due to this fact we are interested to study a theory 
that for low energies, we shall assume reduces to the Maxwell one in the limit of small fields. 

Gravitoelectromagnetic Inflation (GI) was developed very recently with the aim to describe, in an unified manner, 
the inflaton, gravitatory and electromagnetic fields during inflation[l£, 20]. In this formalism all the 4D sources have 
a geometrical origin. This formalism can explain the origin of seed magnetic fields on cosmological scales observed 
today. This proposal was constructed from a 5D vacuum state on a bcd = globally flat metric. As in all 
Space Time Matter (STM) models[2l|, the 4D sources are geometrically induced when we take a foliation on the fifth 
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coordinate which is spacelike and noncompact. However, in the previous works was used the Feynman gauge in order 
to simpHfy the structure of the field equations. 

In this letter we shall use this formalism using Aq = 0. As we shall see, the field equations become coupled, which 
has interesting physical consequences. We shall study the origin and evolution of the seed large-scale electric and 
magnetic fields in a Aq dominated early universe, from a 5D vacuum state, where the expansion of the universe is 
driven by the inflaton field. 



II. VECTOR FIELDS IN 5D VACUUM 



We begin considering a 5D manifold A4 described by a symmetric metric gAB — gsA^ ■ This manifold M is mapped 
by coordinates {a;'^}. 



which, we shall consider as Riemann-flat Rg/jD = 0. To introduce the fields we can define an action in M. 
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R is the 5D scalar curvature. We shall consider these fields as minimally coupled to gravity. In this space the fields 
are free of interactions. The Faraday tensor is antisymmetric Fbc = V^y^c — VcAb- 



A. The 5D Riemann-flat metric with decaying parameter 



In particular, in this letter we are interested to deal with the following Riemann-flat metric [2! 
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where dr^ — dx^Sijdx^ is the euclidean line element in cartesian coordinates and ip is the space-like extra dimension. 
Adopting natural units {h = c = 1) the cosmological parameter A{t) (with A < 0), has units of (length)~^ . The metric 
([3]) is very interesting to study the evolution of the gravitoelectromagnetic (vectorial) field, because is Riemann-flat, 
but has some connections F*^^ ^ 0. This fact is very important when we consider the covariant derivative of . 
The equations of motion for the components of the vectorial field A, are 



d_ 

di 



(v.1) 

d^A, 



A 
— ( 

3 



A 
2A 



dAj _ 
d fdAi 



3 3 oyj \ / 



A 
2A 



dt V dip 
OA, A 

d_ 

dx^ 



2ipe 



2/^ 



}dt 



dAj 



dt 



V'A, 



A 



2 dAi 

tp dip 



0, 



2-tP^ 
3 dx^ 



= 0. 



(4) 
(5) 



(6) 



These are our equations of motion on the metric ([31), once we consider the gauge A" = 0. To solve these equations we 
can begin considering V ■ A = f{t, x, V'), and A4 = ip{t, x, ip). Next, we make Fourier transforma in eqs. (|4]) and (O, 
and use separation of variables for both, fk{t,ip) ~ Fi{t)F2{ip) and (pk{t,ip) ^ a{t)j3{tp) -where we drop de subindex 
k for the transformed variables Fi, F2,a, (3. Working out ^ and ([5]), we arrive to 



^ In our conventions capital Latin indices run from to 4, greek indices run from to 3 and latin indices run from 1 to 3. 



3 



b'd = — Fi, 



a — rd + fc^TTpra 
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where primes and dots denote respectively the derivatives with respect to ^ and t. Furthermore, we define b{t) = 

gl '^'^V^ and T{t) = ^ — \/3A. The constants Ai, A2 come from the separations of variables procedure. To obtain a 
and Fi, we work with the equations ([7]) and ([5]). We introduce a = a and we have 
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where the solutions are given by the primitives 



CT{t)dt, 



a{t) = 

Fl{t) = J CT{t)b{tfdt. 

To solve the equation for F2 we replace (fTO|) and its derivative in ^ 



^P'^F^' + 2ijF^ + A1A2 = 0. 



The solution is 
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i — A1A2. In order to illustrate the formalism we can study an example, which is interesting for the 



with w = 

cosmological expansion of the early universe. 



III. AN EXAMPLE WITH A = Aq: DE SITTER EXPANSION 

We consider the case where A = Aq. When we make the foliation tp ^ ipQ = — F[q^ on the ([3]), this case 

give us a de Sitter inflationary expansion of the universe with tetra-velocities: u" = (1, 0, 0, 0) for a comoving frame. 
Furthermore, the effective 4D line element is 

c/s2 = dt^ - Hq^ e'^"°'dr^. (17) 

In this case equation pT|) yields 

CT + 5Ho& + Hi (/c^e-^^o* _^ g _ ^^^2) CT = 0, (18) 



which has the solution a{t) = e^t^o* [iVi (A;e"^"*) + 7V2>"-ai (^e^^o*)] , with ^ = Y'i + AiA2. To get a{t) 
we have to integrate the Bcssel functions with the exponential, this can be achieved by changing variables to 77 = 
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Hq J dte 2 ^0* = e 2 ^at ^ 'pj^g primitive of a first kind Bessel function gives a Regularized Hypergeometric Function 



(RHF), and for the second kind Bessel function we obtain the solution as a combination of two RHFs 



Hq \ 4, J 

^M22~<^'+^^fcV+§^r(5+^2>.)iF2 



Ho 



5-2/i 
4 

5 + 2n 



; |1 - 



9-2/i 
4 

9 + 2^ 



(19) 



where Mi = iVi - (^^^riTy and A/2 = iV2. To calculate Fi we can define r]' = -\Hq J dter^"°^ = e"^^"* 
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The total solution ipk(t, tp) = [(f'l^^°"^\t) + ak{t)]f3{4>), where we have included the homogeneous solution 



(21) 



that has the typical scale invariant spectrum of a de Sitter model. 

Once obtained the solutions for ipi^:{t,'ip) and fkit^ip), we can try to solve the equations for the potential 3-vector 
Aj{x"^). We take the Fourier transform in the x-space, as before. From eq. ([6]), we obtain the following equation for 
the modes Ck\t,tlj): 

C'U) I TT M) I tt2 ,2 -2Hot A]) j,2 717-2 

4fe +/iotfe +H„k e kk ~V Hq 
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where lC{t, ip) is the source term and 



Mt) = ^[^''"'\t)+ak{t). 



(22) 



(23) 



In this letter we consider that the particle excitations for an observer in ([3|) appear are the Mellin transform in the 
extra coordinate ip. We will see that the extra terms will become massive terms for each m-mode. Thus, this extra 
coordinate formalism, besides it produces couplings between the effective vector and scalar components of the field 
(in a curved background), provide us of a contribution to the mass of vector excitations. The Mellin transform on the 
foliated spacetime (fT7|) . is 
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and the equation ([22]) becomes 

C + HoSj^l + Hi [k'e-'""' - m{m - 1)] = -zfc,H^/C™(t), 



(25) 



where m is a free parameter on the metric (|17p [but not on the 5D Riemann-flat metric ([3])], to be experimentally 
determined by the spectrum of large scale magnetic fields. The total solution to this ordinary differential equation is 
the homogenous part plus the inhomogeneous one 
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The source term, after the Mellin transform, yields 
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Notice that Fi(r) and a(T) are given by hypergeometric functions and we have to integrate them with the Bessel 
functions. This can be done analyticaUy by evaluating the integral for each term of the power series of the RHF. To 
make it, we write explicitly the RHFs in the form 
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The integrals we have to evaluate are of the form 
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where Ki_j-^[x(t)] can be either the first kind or second kind Bessel function. To solve the integral again we repeat the 
procedure of changing variables ??± = — (±Ai + |) / dre^^'^'^^i^^"'^ ~ g-(±A'+i)^for ^ integrals in ([30)1 reduce 

The result for the inhomogencous solution (p7)l . is 
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where //i = —^2 = A* ^ind wi = —W2 = w. 

To study the modes of cosmological interest, we have to work with those modes that during inflation stay outside 
the horizon. This is the limit ke^^°^ ^ 1. Then we truncate the power series to the first term p = q = {). 
We also introduce a constraint in the possible values of m. If we look at equation ((25)) . we see that, in order to 
identify this term with ordinary matter, the quantity m{'m — 1) > in eqs. ([5^ and This means that to < or 
TO > 1. Then, to keep positive parameters for the Hankel functions in (j26p . we shall work with the negative values 
of TO. Notice that and w are not independent parameters. They are related by the equation, + uP' 
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we require that fi and w be real and positive, then we obtain that they are restricted to the interval 
restricts the parameter space of Ai, A2. 
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Considering that the magnetic fields produced during this epoch are scale invariant in the cosmological level, we 
can define a value for the parameter m 



1 f6kH(t) Ji. 
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The contribution to the magnetic fields is exclusive from the homogeneous solution in (|26|) . this should be clear from 
^ because the source terms are all gradients. The k-power that come from the homogeneous solution of the modes is 
m — ^. Then, to obtain a nearly scale invariant spectrum of the magnetic field we need that 5 — 1 + 2m ~ 0, yielding 
TO ~ — 2. In the cosmological limit the inhomogeneous equation can be reduced to 
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where e is a parameter that takes into account small deviations from scale invariance. In the last expression should be 
noted that we have drop the terms involving positive /x-powers of the physical wavenumber fce"^"* (this because of the 
same idea we stayed with p = q = 0). The Bessel functions have been approximated to their asymptotic expressions 
in the infrared limit. Note that the term involving e decays very strongly due to the factor x{t)^. 

A. Quantization and normalization of homogeneous solutions of the vector and scalar fields 

We consider the usual commutation relations for the fields and their conjugate momenta on the effective 4D metric 
(fT7)) . For the scalar field we obtain 



[^™(t, f), (^,„(t, f )] = ^i^3e-=^^«*(5(3) (f - f ), 
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where fmit, x) = ip{t, x)Prm E^nd is the Mcllin transform of I3{tp). From these relations wc derive the normalization 
condition for the modes 



(41) 



where the Bunch-Davies vacuum is given by Ai = and then, we obtain Ao — i ^"V^ . We repeat the same for the 
vector solution, the commutation relation is in this case 

A™j(i,f),i™-,fe(i,f )] = iSjkSmrn' Hoe-''<'*S^''\x - x'), (42) 

where Am,j{t, x) is the Mellin transform of the Aj field. We impose the Bunch-Davies vacuum to the modes choosing 
Di = and D2 = in ((26|) . Then, the Fourier-Mellin modes comply 
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We have noted that the magnetic fields depend exclusively from the homogeneous solution, then we can compute their 
quadratic amplitude from (|38p 
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The horizon wavenumber is found to be kH{t) = fn{m — l)e^°^. This means that the physical magnetic field is 

constant for any m 
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(277)3/22™ ° 



0ji + m(m-l) 



2+m 



(46) 



which is divergent for m = —2. If we require that these fields are nearly invariant in cosmological scales, then 
m = — 2 + e, e being a small parameter, this yields for the amplitude of the physical magnetic field 
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On the other hand the electric fields are affected by the extra terms. Then, at the end of the inflationary epoch there 
will be a contribution from the inhomogeneous solution to the spectrum and amplitude of the electric field. This field 
has components E" = F°''^V0, where are the components of the observer velocity. For a physical observer we have 
Ephys = a~^{t)dtA. Therefore 
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Prom the homogeneous solution product we obtain the amplitude 
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Hence, the homogeneous amplitude of the electric field, is 
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The relation between the energy density of the electric and magnetic fields, in the physical frame, is 
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For a nearly scale invariant magnetic field, we obtain 
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This means that the energy of the electric field is dominant during exponential infiation on cosmological scales. For 
a nearly scale invariant electric field we obtain m = — 3/2 + e, and then 
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Thus, the magnetic field would be constant but bigger on smaller (astrophysical) scales, while the electric field would 
still dominate in larger scales 
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B. Spectrum of the inhomogeneous solutions of the electric field and the scalar 



The amplitude and spectrum of the electric field have terms that involve the inhomogeneous solution, with double 
infinitum power series, coming from an integral of the hypergeometric functions. For convenience we only keep the 
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first term {p = q = 0), because it is of cosmological relevance. The contribution of terms with p and q can be 
neglected on cosmological scales. The three remaining contributions are 

2, _ Hje-^""' f dk 3■[h]U)■[^„h]Uh 

\^phys I — / 'JT ^k,m ^k,m ' 

/p[2] 2, _ H^e-^"»* f ji-linh]U)^lh]U)* 
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\^phys I ^ 27r2 / k '"^ '™ ' 

In order to simplify the notation, we write (|39p in the compact form 

t:^{t) = -i|ie-^°* ^ Di":),^, (58) 

n=l,2 

where we used kj = kej and the coefficients ©^"^ have units of Hq. This yields the respective power spectrums, for 
(USD, dSSl) and 

/f (60) 

where we have considered only terms with p = g = in eqs. (|32|) and ()33p . It is important to notice that the spherical 
symmetry is broken because ej is an unitary vector. 

On the other hand, the power spectrum for the homogeneous part of the scalar field is scale invariant for the solution 
([23| . while that for the other terms in ([23| we obtain two different spectrums to zero order in the hypergeometric 
function (gH) 

- ^/ fk^ ^^t"'Ht)aUt) - e-(5/-M) J (61) 

= i / f fc' mi'"""^ (tMt) - I f (62) 

('/'[^l') = i / f Mt)am - e-2(^/2-'') I ^ (63) 

Notice that these inhomogeneous terms are exponentially damped. The parameter fi can be fixed, with the smaller 
index that decays weaker, so as to yield the experimental data for the scalar spectral index [23[, rig — 0.958, then we 
write 3/2 - ^ = e', with e' = Us - 1 0.042. 

IV. FINAL COMMENTS 

In this letter we have studied the primordial spectrum of electromagnetic fields using GI. Starting from a gauge 
with ^0 = 0, we have obtained some interesting properties. In the example here studied the spectrum of large-scale 
magnetic fields is nearly scale-invariant for m ~ —2. The amplitude for the strength of comoving magnetic fields is 
dramatically increasing, but they are frozen in physical coordinates. The important result here obtained is that the 
modes of Aj are affected by a source, which is originated in the modes of the inflaton field, so that the spectrum of 
the large-scale electric field during infiation depends of the modes of the infiaton field. These modes can be considered 
as massive photons which are gauge-invariant in a 5D sense, but once the foliation dip = is done (which implies the 
choice of a relativistic system), these photons acquire mass because they live in an effective 4D curved spacetime. In 
this sense the choice of the relativistic system acts as an effective Higgs's mechanism. 

But the more interesting result relies in that the spectrum of the infiaton field depends of the modes of Ai, because 
they are coupled to the modes of Ai. Of course, this scale invariance is significatively affected on shorten scales, so 
that it is nearly scale invariant on very large scales. This result disagrees with standard 4D versions of inflation, but 
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it agree very much with experience, because it is very known that for shorten scales the mass spectrum of matter has 
a positive index with a scale dependent power. 

The effects of a conducting plasma in the early inflationary universe are negligible. During inflation conformal 
invariance is broken and the strength of comoving magnetic fields increases dramatically as a} , until values of the 
order of Bcom — 10^^^ Gauss after 63 e— folds, so that the flatness problem is resolved in the model. After inflation, 
the universe enters in the so-called reheating phase, during which the energy of the inflaton is converted into ordinary 
matter. In this epoch, the conductivity of the universe is of the order of <7c ^ T ^ H (with a background 
temperature T <C Mp). magnetic fields evolves adiabatically from the end of inflation until today, due to the high 
electrical conductivity of the cosmic plasma. In this epoch the universe is thermalized, so that the comoving magnetic 
field decreases with the expansion to take actual values of the order of 10~^ Gauss[l^, Notice that the results 
here obtained depends on the gauge Aq = 0. It is well known that any viable mechanism to generate seed magnetic 
fields during inflation must repose on the breaking of conformal invariance of standard electrodynamics. Otherwise, 
the produced fields are vanishingiy small. Notice that the approach here worded is not conformally invariant on the 
effective 4D metric ([TTl) . The origin of this rupture is in the fact that some connections are non-zero on the 
5D Riemann flat metric ([3|). This is the reason by which bosons are massive on the effective 4D spacetime ([T7| on 
which move the observers. Concerning electric fields, there is a damping of the longitudinal component of the field 
strength, corresponding to the gradual neutralization of charged particles in the primordial plasmafHI, in the first 
stages of reheating. Finally, in our model, inflation occurs at a very low scale with H ~ 10~^ Mp and with the inflaton 
field taking values much below of the Planckian scale: {(p) ~ 10"^^ Mp26]. In this sense our model evolves on scales 
similar to the MSSM inflationary model(27j. where fine tunning and slow rolling problems joined with reheating were 
considered and the inflaton field couplings to Standard Model physics is explained from flrst principles. In our case 
the couplings between the fields Ac is explained from the induced curvature of the metric P7|). The problem of 
back-reactionl28i should be considered in future works. 
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